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Abstract 


Three  generations  of  computers  have  elapsed  since  the  first  satisfactory 
method  for  solving  transportation  and  transhipment  problems  was  devised. 
During  this  time  many  computational  advances  have  taken  place  in  developing 
computer  codes  to  solve  these  preplans.  The  primary  purpose  of  this  paper  is 
to  sunnarize  these  events  and  to  do  sane  crystal  ball  gazing  to  provide  what 
we  believe  to  be  "best  estimates"  of  future  trends. 
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PAST 


Approximately  200  vears  have  elapsed  since  the  French  Academy  of 
Sciences  posed  the  civil  engineering  problem  of  'cutting  and  filling." 

Their  formulation  of  the  problem  was  not  the  same  as  the  transportation 
problem  as  we  know  it  today,  but  was  the  equivalent  of  a  transportation 
problem  in  continuous  form.  The  current  formulation  of  a  transportation 
problem  was  due  to  Kantorovich  (1939) ,  Hitchcock  (1941) ,  and  Koopmans 
(1947). 

Kantorovich  showed  in  1939  that  a  class  of  problems  closely  related 
to  the  transportation  problem  had  a  remarkable  variety  of  applications. 

These  were  concerned  typically  with  the  allotment  of  tasks  to  machines  whose 
costs  and  rates  of  production  vary  by  task  and  machine  type.  Kantorovich 
gave  a  useful  lait  inoanpiete  algorithm  for  solving  such  problems.  Later, 
in  1942,  he  studied  both  discrete  and  continuous  versions  of  this  problan 
and  in  1948,  along  with  Gavurin  wrote  an  applicational  study  on  the  capa¬ 
citated  transportation  problan. 

Hitchcock  developed  an  inconplete  algorithm  in  1941,  which  exploited 
special  properties  of  the  transportation  problem  to  find  starting  solutions, 
/oopmans  (1947)  independently  arrived  at  the  same  problem  in  connection 
with  his  work  as  a  member  of  the  Combines  Shipping  Adjustment  Board.  He 
and  Reiter  discussed  the  problem  from  an  economic  efficiency  analysis 
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viewpoint  and  pointed  out  the  antilogy  between  it  and  the  classical  Maxwell- 

J 

Kirchhoff  electrical  network  problem.  Because  of  their  work,  the  problem 
is  often  referred  as  the  Hitchcock-KQopmans  Transportation  Problem. 

The  first  generally  satisfactory  method  for  solving  the  general  class 
of  transportation  models  was  due  to  G.B.  Dantzig  in  1949.  ihis  method  is 


sometimes  called  the  Row-Column  Sun  Method  [  4  ]  or  the  MODI  method  [7] . 
Chames  and  Cooper  (1954)  later  wrote  an  explanation  (dubbed  the  Stepping- 
Stone  Method)  of  the  simplex  steps  involved  in  the  Row-Column  Sum  Method.  With 
the  advent  of  a  method  for  solving  the  transportation  problem  came  numerous 
methods  for  securing  starting  bases.  Two  of  the  methods  ccmnonly  referenced 
are  the  Northwest-Comer  Rule  [  7  ]  and  the  Vogel  Approximation  Method  [  31  ] 
(often  referred  to  as  VAM) .  Of  all  the  start  methods  developed,  VAM  became 
the  one  most  used  for  hand  calculations  due  to  the  excellent  start  it  pro¬ 
vides.  Thus  in  the  folklore  VAM  is  considered  the  best  procedure  for  both 
ocrputer  and  hand  calculations. 

After  the  development  of  the  Rcw-Colurnn  Sum  Method,  the  transportation 
model,  with  integer  parameters,  rapidly  became  the  chief  "solvable"  integer 
linear  programing  problem  due  to  the  integer  extreme  point  property.  Also 
a  survey  by  L.W.  smith,  Jr.  in  1956  indicated  that  at  least  half  of  the  linear 
programming  applications  used  this  model.  Some  reasons  for  the  surprising 
concentration  on  problems  of  this  kind,  particularly  in  applications,  are: 

(1)  Answers  to  "large"  problems  can  be  easily  computed  by  hand, 

which  is  an  impossible  task  for  general  linear  programming  problems  of  srnilar 
dimensions.  Also,  integer  solutions  were  immediately  attainable. 

(2)  It  is  possible  to  approximate  many  linear  programs  by  transportation 


problems . 


(3)  A  number  of  seemingly  unrelated  linear  programs  have  been  found  to 
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be  equivalent  to  transportation  problems. 


(4)  Conputer  codes  were  developed  as  early  as  1952  for  solving  trans¬ 


portation  problems. 


(5)  Business  executives  can  understand  the  transportation  model,  leading 


to  increased  demand  for  its  applications  in  practical  settings. 


Subsequent  to  the  development  of  the  simplex  based  Row-Column  Sum 


Method,  Ebrd  and  Fulkerson  (1955) ,  developed  a  primal-dual  method  for  solving 


transportation  problems.  Somewhat  earlier  Gleyzal  [18]  developed  a 


method  similar  to  the  primal-dual  method. 


It  is  interesting  to  note  that  Dantzig  and  Ford  and  Fulkerson  concluded 


on  the  basis  of  hand  calculations  that  the  primal-dual  method  was  superior  in 


efficiency  to  the  Row-Column  Sun  Method.  This  conclusion  was  also  sup¬ 


ported  by  Flood  (1961)  on  ocnputer  codes.  Consequently,  this  conclusion 


became  part  of  the  folklore. 


The  first  computer  code  for  solving  transportation  problems  was  based 


on  the  Row-Colunn  Sun  Method  and  in  terms  of  current  jargon  is  called  a 


primal  transportation  code.  About  1952  such  a  oode  was  developed  by  the 


George  Washington  university  Logistics  Research  Project  in  conjunction  with 


the  Confutation  Laboratory  of  the  National  Bureau  of  Standards  [33  ].  This 


code,  which  was  designed  for  use  on  the  Bureau  of  Standards  Eastern  Auto¬ 


matic  Conputer,  was  further  improved  try  the  NBS  Confutation  Laboratory.  The 


code  was  capable  of  solving  problems  with  at  most  600  nodes  with  a  pivot  time 


of  3  plus  minutes  pei  pivot.  Current  pivot  times  for  problems  of  this  size 


are  6-10  milliseconds  on  the  CDC  6600,  UNIVAC  1108,  and  IBM  360/65  computers. 


The  inoore-out-of-core  primal  code  by  Dennis  (1958)  is  one  of  the  first 


codes  to  be  described  in  detail  in  the  literature.  Dennis’  paper  is  also  one 
of  the  first  to  study  different  criteria  for  selecting  pivot  elements.  Un- 
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fortunately,  his  study  principally  involved  onl^  one  problem  of  size  30 
origins  by  260  destinations.  The  best  solution  time  was  9.6  minutes  on  the 
Whirlwind  ccrputer. 

Another  primal  code  was  developed  in  the  Tate  fifties  by  Glicksman, 

Johnson  and  Eselson  (1959) .  This  code  was  developed  for  the  UNIVAC  I  for 
solving  "thin  rectangular"  problems  and  was  also  an  incore  -  out-of-oore  code. 

The  code  solved  a  15  origin  by  488  destination  problem  in  24  minutes.  This 
is  approximately  200  times  slower  than  current  incore  codes. 

Also,  during  the  late  fifties  a  number  of  transportation  codes  were  de¬ 
veloped  using  Kuhn's  Hungarian  Metliod,  implemented  primarily  on  IBM  computers. 
These  codes  include  the  one  due  to  Flood  (1961)  using  his  proof  of  the 
Konig-Egervary  Theorem.  Codes  based  on  Ford  and  Fulkerson's  primal-dual 
algorithm  were  also  beginning  to  be  implemented,  such  as  the  IB-TFL  code 
(1958) .  The  code  of  Flood  and  the  IB-TFL  code  were  compared  on  a  problem 
with  29  origins  and  116  destinations  on  an  IBM  704.  Their  times  were  3.03 
and  3.17  minutes,  respectively.  Current  solution  time  would  probably  be  2  seconds. 

Following  these  developments,  there  was  a  histus  of  half  a  dozen  years 
during  which  little  was  visibly  aocatplished  in  the  development  of  improved 
solution  methods  or  computer  implementations.  From  an  algorithmic  stand¬ 
point,  it  was  widely  believed  that  no  significant  refinements  remained  to  be 
discovered.  In  retrospect,  this  attitude  seems  surprising,  particularly  in 
view  of  the  paucity  of  experimentation  to  determine  the  computational  strengths 
and  weaknesses  of  alternative  approaches.  Then,  in  the  later  sixties  and 
more  particularly  in  the  early  seventies,  a  new  surge  of  interest  in  network 
methods  and  applications  came  about,  leading  to  a  number  of  surprises  for 
those  steeped  in  the  notions  of  a  decade  earlier. 

It  is  to  these  more  recent  developments  t  at  we  now  turn. 
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PRESENT 


As  already  intimated,  the  early  special  purpose  primal  and  primal-dual 
codes  were  capable  of  solving  only  snail  problems,  were  quite  slow,  and 
were  not  extensively  tested.  Begin,  ling  with  the  latter  half  of  the  sixties, 
several  codes  have  been  jointly  developed  by  mathematical  programmers  and 
systems  analysts  who  have  performed  extensive  experimentation  on  various 
algorithmic  rules.  In  particular,  recent  code  development  and/or  compu¬ 
tational  studies  have  been  performed  by  Barr,  Glover,  and  Klingman  (1972) 
on  an  out-of -kilter  network  code,  Bennington  (1971-72)  on  a  non-extreme 
point  network  code,  Boting  (1966)  on  an  out-of-kilter  code,  Clasen  (1966) 
on  an  out-of-kilter  code.  Control  Data  Corporation  (1970)  on  an  out-of¬ 
kilter  code,  Glover,  Kamey ,  Klingman,  and  Napier  (1970-1972)  on  a  primal 
transportation  code,  Glover,  Kamey,  and  Klingman  (1970-1972)  on  a  dual 
transportation  code,  Glover,  Kamey,  and  Klingman  (1973)  on  a  primal  tranship¬ 
ment  code,  Klingman  Napier,  and  Ross  (1973)  on  parameter  sensitivity 
analysis,  Klingman,  Napier,  and  Stutz  (1973)  on  a  network  generator  code,  Lee 
(1968)  on  various  solution  procedures,  Srinivasan  and  Thompson  (1970-1972) 
on  a  primal  transportation  code,  Texas  Water  Development  Board  (1961)  on  an 
out-of-kilter  code,  and  UNIVAC  (1973)  on  an  out-of-kilter  code.  It  should  be 
noted  that  all  of  these  codes  are  in-oore  codes  coded  in  FORTRAN. 

The  primal  transportation  and  the  primal  transhipment  codes  by  Glover, 
Kamey,  and  Klingman  [13,  153  utilize  the  augmented  predecessor 
index  r.ethod  [  14  ],  which  elaborates  on  Johnson's  "triple-label  method"  [20  3 
by  providing  an  efficient  method  for  characterizing  successive  basis  trees 
with  minimal  relabeling.  The  augmented  predecessor  list  structure  has  been 
a  major  contributor  to  the  improvements  in  the  ocmputational  efficiency  of 
solution  algorithms.  With  its  use,  the  primal  transportation  cod  >  by  Glover, 
Kamey,  and  Klingman  [  15  ]  executes  a  pivot  on  a  600  node  problan  in  6 
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milliseconds  compared  with  the  early  breakthrough  (1952)  of  reducing  the 
time  per  pivot  to  3  minutes.  While  this  reduction  is  largely  due  to  improve¬ 
ments  in  computers  and  the  in-core  nature  of  these  codes,  this  is  not  the 
whole  reason.  Fbr  instance,  the  first  accelerated  primal  transportation 
code  developed  by  Srinivasan  and  Thompson  [  34  ]  employed  a  list  structure 
for  preceeding  in  a  forward  direction  through  a  spanning  tree  similar  to 
Dennis'  procedure.  Upon  comparing  solution  times  of  the  Srinivasan  and 
Thompson  code  with  the  Glover,  Kamey,  and  Klingman  oode  on  the  same  problems 
and  machine,  the  efficiencies  of  the  augmented  index  structure  became  ap¬ 
parent.  Srinivasan  and  Thorp  son  recoded  using  the  augmented  list  structure 
and  cut  their  solution  tunes  by  more  than  half. 

The  oode  development  and  comparison  of  Srinivasan  and  Thompson  [34  ] 
provides  an  important  ocrputational  analysis  of  several  primal  start  proce¬ 
dures  and  pivot  criteria.  The  purpose  of  this  study  was  to  determine  a 
design  for  an  in-core  uncapacitated  primal  transportation  and  assignment  oode 
which  optimally  combines  start  procedures  and  pivot  criteria  for  maximum 
solution  efficiency-  The  study  disclosed  that  the  best  start  method  is  the 
"modified  row  minimum  start"  procedure  and  the  best  pivot  selection  criterion 
is  the  "row  most  negative  rule."  This  pivot  rule  was  also  found  to  be  best 
by  Dennis  [8  ].  Maximum  problem  size  solved  was  350  nodes  (origins  plus 

destinations).  This  node  limitation  is  due  to  the  fact  that  the  code  was 
all  in-core  and  stored  a  complete  cost  matrix.  The  average  solution  time  on 
175  origin  by  175  destination  transportation  and  assignment  problems  was  7.8 
seconds . 

The  code  development  and  comparison  by  Glover,  Kamey,  Klingman,  and 
Napier  (1970-1972)  performed  similar  analyses  on  a  broad  profile  of  dense 
and  nondense  problems.  The  underlying  algorithm  was  specially  designed  for 
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solving  both  capacitated  and  uncapacitated  problems  with  nondense  cost 
matrices  (i.e. ,  transportation  problems  where  sane  cells  may  not  be  allow¬ 
able)  .  This  study  also  found  the  modified  row  minimum  start  and  row  most 
negative  pivot  rule  to  be  best,  thus  casting  doubt  on  the  folklore  of  the 
superiority  of  VAM  starts.  In  addition,  the  study  compared  100  problems  on 
this  primal  code  to  several  other  codes  including  Clasen's  SHARE  code  (1966), 
the  Glover,  Kamey,  and  Klingman  dual  code  (1970),  and  the  state-of-the- 
art  linear  programming  code  OPHELIE/LP.  This  caparison  revealed  that  the 
specially  designed  primal  code  was  at  least  eight  times  faster  than  the 
SHARE  and  dual  codes,  and  150  times  faster  tlian  OPHELIE/LP.  Thus  the  old 
folklore  about  the  superiority  of  out-of-kilter  methods,  and  a  new  folklore 
among  computer  service  divisions  about  equivalence  of  general  purpose  and 
special  purpose  solution  codes  for  transportation  and  transhipment  problems 
were  upended. 

The  largest  problems  solved  by  the  Glover,  Kamey,  Klingman  and  Napier 
study  were  1000  origin  by  1000  destination  problems  with  an  average  solution 
time  of  17  seconds.  This  study  also  tested  the  primal  code  on  four  computers, 
IBM  360/65,  UNIVAC  1108,  CDC  6400,  and  CDC  6600  in  order  to  provide  insights 
into  conclusions  based  on  oarparing  times  on  different  machines  and  com¬ 
pilers.  It  was  discovered  that  standard  guidelines  concerning  the  relative 
efficiencies  of  different  computers  were  completely  misleading,  since  the 
primal  code  ran  only  10%-12%  faster  on  the  CDC  6600  than  on  the  UNIVAC  1108 
and  the  IBM  360/65  differing  substantially  from  the  estimates  one  would  ob¬ 
tain  by  comparing  instruction  execution  times  of  the  machines. 

Motivated  by  the  fact  that  out-of-kilter  codes  were  found  to  be  substan¬ 
tially  slower  than  the  special  primal  code  ,  Barr,  Glover,  and  Klingman  (1972) 
developed  an  improved  version  of  the  out-of-kilter  method  which  was  sub- 


sequently  coded.  This  code  was  found  to  be  only  40%  slower  (on  the  same 
problems  and  machine)  than  the  primal  transportation  code  of  Glover,  Karoey, 
and  Klingman  on  transportation  problems.  This  code  was  also  ccrrpared 
against  Clasen's  SHARE  code,  Boeing's  code,  and  the  Texas  Water  Development 
Board  code  and  found  to  be  at  least  six  times  faster  than  the  best  of  these 
(which  differed  from  problem  to  problem) .  The  study  also  examined  a  total  of 
215  capacitated  and  uncapacitated  transhipment  problems  demonstrating  the 
superiority  of  the  improved  version  of  the  out-of -kilter  code  over  the  other 
out-of-kilter  codes  in  all  cases.  The  largest  oroblems  solved  were  1500 
node  transhipmant  problems.  The  mean  solution  time  was  34  seconds. 

Still  more  recently,  Glover,  Kamey,  and  Klingman  (1973)  have  developed 
a  general  primal  transhipment  code.  Computational  conparison  of  thi^  code  with 
the  out-of-kilter  code  by  Barr,  Glover,  and  Klingman  reveals  that  the  primal 
code  is  30%  faster  on  transhipment  problems.  This  is  rather  startling  since  the 
Barr,  Glover,  and  Klingman  code  is  probably  the  fastest  out-of-kilter  code 
in  the  world  and  conventional  wisdom  has  it  that  labeling  techniques  are 
inherently  more  efficient  than  simplex  techniques.  The  primal  transhipment  code 
was  also  tested  against  an  out-of-kilter  code  by  Bennington  (1971-72)  and 
found  to  be  eight  times  faster.  This  computational  study  also  shows  the  super¬ 
iority  of  the  new  primal  transhipment  code  in  terms  of  central  memory  require¬ 
ments  for  storing  network  data.  Specifically,  the  out-of-kilte.’*  codes  dis¬ 
cussed  earlier  require  6-10  arc- length  arrays  and  4-7  node- length  arrays 
as  compared  to  3  cure- length  arrays  and  8  node- length  arrays  for  the  primal 
oode.  The  substantially  increased  problem  size  that  can  be  acocitmodated  by 
the  new  primal  code  i3  illustrated  in  the  study  by  the  solution  of  an  8000 
node  problem. 


In  addition  to  recent  code  development  efforts  and  ocnputational  com¬ 
parisons,  two  computational  studies  on  the  effects  of  parameter  values  have 


!) 


been  conducted  on  transportation  problems.  The  first  study,  by  Srinivasan  and 
Thompson  (1970-71) ,  examined  the  effect  on  solution  time  of  rect angularity , 
the  mmber  of  significant  digits  in  uniformly  distributed  cost  coefficients, 
and  the  mxnber  of  significant  digits  in  the  supply  and  demand  values.  The 
study  analyzed  600  transportation  problems  with  a  maximum  size  of  350  nodes. 

The  main  finding  is  that  solution  time  is  not  particularly  sensitive  to  the 
supply  and  demand  ranges ,  but  is  quite  dependent  on  cost  ranges.  The  second 
study,  by  Klingman,  Napier,  and  Ross  (1973) ,  performed  a  detailed  examination 
of  the  effects  of  problem  dimensionality  on  solution  times.  The  study  included 
over  1000  randomly  generated  problems  with  185  different  combinations  of  number 
of  origins,  number  of  destinations  and  number  of  variables  (not  all  cells  being 
considered  admissible) .  Every  problem  was  solved  using  three  starting  pro¬ 
cedures.  Over  10,000  pieces  of  data  were  analyzed,  providing  numerous  insights 
into  the  computational  effects  of  the  number  of  constraints,  the  degree  of 
" rectangular! ty , "  and  the  number  of  variables.  An  unexpected  finding  of  the 
study  was  that  no  single  starting  procedure  dominates  another;  rather,  the 
efficiency  of  a  starting  procedure  varies  with  problem  structure. 

These  studies  on  the  effect  of  parameter  values  and  code  comparisons 
repeatedly  reinforced  the  conclusion  that  in  order  for  researchers  to  compare 
their  codes  in  a  meaningful  way  it  is  necessary  that  they  use  exactly  the 
same  problems.  This  is  due  to  the  fact  that,  even  if  two  randomly  generated 
problans  have  the  same  parameter  values,  a  generator  inherently  builds  structure 
into  the  problems,  particularly  transhipment  and  transportation  problems  in 
which  seme  arcs  do  not  exist.  This  point  is  underscored  in  the  computational 
study  by  Barr,  Glover,  and  Klingman  and  has  been  further  demonstrated  by  the 
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comparison  of  the  codes  due  to  Bennington;  Glover,  Klingnan,  and  Kamey;  and 
Srinivasan  and  Thcnpson  which  yielded  unexpected  results  when  tested  on  the 
same  problems  and  the  same  computer. 

To  enable  researchers  to  ntaningfully  oonpare  their  solution  codes, 

Klingman,  Napier  and  Stutz  [24]  developed  a  code  which  generates  assignment  prob¬ 
lems,  and  capacitated  and  uncapacitated  transportation  and  transhipment  problems. 
In  addition  to  producing  structurally  different  classes  of  transhipment  problems, 
the  code  permits  the  user  to  vary  structural  characteristics  within  a  class. 

By  means  of  this  code,  researchers  can  generate  identical  transhipment  problems 
(independent  of  the  computer).  Advantages  of  the  transhipment  generator,  which 
is  available  with  documentation  to  researchers  for  a  nominal  handling  charge, 
are  its  ease  of  use  (requiring  only  two  data  cards  per  transhipment  problem)  and 
the  standardization  of  its  output  (generating  problems  for  use  by  other  codes 
in  SHARE  input  format).  In  addition,  the  latter  part  of  the  documentation 
provides  the  user  with  the  data  on  40  assignment,  transportation,  and  tran¬ 
shipment  problems  varying  in  size  from  200  nodes  to  8,000  nodes  and  from 
1,300  arcs  to  35,000  arcs.  The  objective  function  value  and  solution  time  for 
these  problems  are  also  provided  for  Clasen's  SHARE  code,  Boeing's  code,  Barr, 
Glover  and  Klingman  out-of-kilter  code,  and  the  Glover,  Kamey  and  Klingman 
primal  transhipment  code. 

FUTURE 


In  spite  of  t  ie  major  recent  gains  in  the  development  and  testing  of 
network  codes ,  significant  avenues  remain  to  be  explored.  In  particular, 
all  the  codes  currently  in  vogue  are  in-core  codes,  all  are  coded  in  FORTRAN, 
most  have  not  fully  exploited  problem  size  capability  of  third  generation 
carputers ,  and  all  of  them  can  only  solve  pure  network  problems  as  opposed 
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to  generalized  or  weighted  network  problems. 

In  the  neai  future,  we  will  undoubtedly  see  large  scale  in-core-out-of 
core  network  codes  developed  which  are  capable  of  solving  network  problems  of 
almost  unlimited  size.  In  this  regard,  there  is  an  in-core-out-of-oore  primal 
transhipment  code  under  development  by  Karney  and  Klingman  which  has  not  been 
fully  streamlined  but  which  has  solved  a  problem  for  the  Naval  Personnel  Labora¬ 
tory  in  San  Diego  with  2400  nodes  and  450,000  arcs  on  the  I  EM  360/65  and  CDC 
6600  in  66  and  65  minutes,  respectively.  Also,  there  is  a  primal  generalized 
network  code  being  developed  by  Glover,  Klingman,  and  Stutz  utilizing  ex¬ 
tensions  of  the  augmented  predecessor  method  [17]  that  has  proved  successful 
for  the  pure  network  codes.  In  addition,  we  shall  probably  see  codes  developed 
in  other  languages  (e.g. ,  ALGOL  and  APL)  in  order  to  rigorously  determine 
which  language  is  best  in  network  applications.  Assembly  language  codes  are 
also  quite  conceivably  in  the  offing. 

While  the  transhipment  generator  [24]  is  a  start  towards  helping  bench¬ 
mark  these  codes,  we  believe  that  a  bureau  needs  to  be  established  to  enable 
standardized  comparisons.  We  have  been  informed  by  Richard  Jackson  that  the 
National  Bureau  of  Standards  in  conjunction  with  various  researchers,  is  con¬ 
sidering  the  feasibility  of  establishing  a  service  facility  to  accomplish  this. 
The  benefit  of  establishing  such  a  service  is  apparent;  however,  numerous  prob¬ 
lems  must  be  overcome.  For  example,  Input/Output  formats  of  codes  and  the  timing 
of  codes  must  be  standardized.  While  such  natters  may  appear  to  be  quite  simple, 
they  are  not.  To  illustrate,  every  code  benchmarked  by  the  authors,  e.g.,  those 
due  to  Clasen,  Boeing,  the  Texas  Water  Development  Board,  General  Motors ,  Ben¬ 
nington,  and  Srinivasan  and  Thompson,  used  a  different  input  format,  and 
considerable  effort  was  required  to  accommodate  these  differences.  Also, 
while  a  valid  criterion  for  the  timing  of  in-oore  codes  is  quite  easy  to 
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define  (namely  centre. processing  tine  exclusive  of  data  input  and  output) , 
fcstabliuhing  an  accepi  able  treasure  for  the  timing  of  in- core -out-of -core 
codes  in  a  multi-proo'issing  environment  is  far  more  complex.  In  any  case, 
we  belxeve  sate  criteria  need  to  be  developed. 

Other  short  ranp  -*  future  developments  which  we  foresee  include: 

a)  development  *;f  network  aonputer  systems  similar  to  general  linear 
programming  systems.  Tliese  systems  will  include  such  things  as  a  cartmand 
language  (which  allows  the  user  to  add,  delete,  and  modify  arcs) ,  matrix 
generators,  report  generators,  user  subroutine  control  of  the  system  or  any 
oonponent,  and  interactive  coupling  with  data  base  management  info  rmation 
systems.  The  Control  Data  Corporation  NLT^LCW  [  29  ]  system  and  the  UNIVAC 
UKILT-1100  [35  ]  system  are  forerunners  of  such  systems. 

b)  establishment  of  numerous  special  purpose  integer  programming  codes 
using  efficient  network  codes  as  the  main  computational  vehicle,  e.g. ,  plant 
location  codes,  fixed  charge  network  codes,  integer  generalized  network  codes, 
constrained  network  codes,  multi-carmodity  (integer)  network  codes,  constrained 
generalized  network  codes,  and  multi-oatmodity  (integer)  generalized  network 
codes. 

looking  farther  into  the  future,  wo  anticipate  the  following  as  possible 
developments: 

a)  an  efficient  graph  compute:  language  which  allows  a  user  to  write 
special  purpose  network  codes  in  half  a  day.  A  forerunner  is  the  GROPE 
language  at  the  University  of  Texas. 

b)  network  and  related  optimization  codes  which  modify  themselves- 
for  example,  computer  network  codes  which  "learn"  how  to  efficiently  solve 
particular  types  of  problems  through  experience  in  solving  than.  (Preliminary 
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investigations  of  this  type  have  now  been  going  on  for  more  than  a  decade. ) 

c)  multi-page  linear  prograrttning  codes  which  use  special  purpose  codes 
(e.g. ,  network  codes)  to  solve  pages  (corrponents)  which  have  a  special 
structure. 
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